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ABSTRACT 


Àn evaluation is performed of a radiation model for the 
Naval Postgraduate School primitive-equation numerical wea- 
ther prediction system. The model employs empirical expres- 
552 .55cor atmospheric absorptivity, scattering-reflectivity, 
cloud-reflectivity and earth-surface reflectivity to compute 
solar insolation absorbed at earth and in the key atmos- 
pheric layers. The terrestrial cooling effect at earth and 
in these same key atmospheric layers is formulated using re- 
cent empirically-derived emissivities for the effects of 


- 


both vater vapor and Co». 
Mean seasonal Nu spheric soundings for the Northern 
Hemisphere are utilized for testing the model. In addition, 
application of atmospheric boundary-layer modeling permits 
determination of the surface-layer turbulent transports (in 
mee vertical) of sensible and latent heat at the earth's 
57:66 To evaluate the validity of the radiation model, 


heat budgets are compiled for the earth-atmosphere system, 


the atmospheric column and the earth's surface. 
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Pa ETC 


The present study was to evaluate a radiation package 
for the Naval Postgraduate School primitive equation pre- 
diction model. The radiation package is WE sett an out- 
growth of the model used in the UCLA general circulation 
model (Gates et al. 1971) with modifications suggested 
from the NCAR model (Oliger et al. 1970) and from the Geo- 
physical Fluid Dynamical Laboratory (Manabe and Strickler, 
1964). 

The model consists of five G-levels (G = P/T where T 
is the surface pressure) for the purposes of data input 
as Bn in Figure 1. Cloud cover was treated by the use 
of a single cloud layer extending from o = 0.8 to 0.4. 

In the computations for both pressure-scaled water-vapor 
.and co, absorber masses, a more realistic method was in- 
troduced to extrapolate accurate water-vapor values above 
levels normally available from radiosonde data (P = 300 mb). 
A system for the integrations of net terrestrial flux was 
introduced by using empirically derived absorptivities 
(Sasamori, 1968) based on the radiation chart of Yamamoto 


(1952), for water-vapor and CO Evaporation and sensible 


2” 
heat fluxes were derived from a formulation for surface 
layer heat transfer by Langlois and Kwok (1969). 


Northern hemisphere climatological atmospheric data 


from London (1957) were utilized to evaluate the model. 





Pressure, temperature, and relative humidity data were 
divided into the four seasons and available over ten de- 
gree latitude bands from the equator poleward. 

The development of the terrestrial fluxes is given 
in Section III and the solar radiation model is given in 
Section IV. The results of the ти long and short 
wave fluxes were then compared to those of Budyko (1955), 
London (1957) and Katayama (1967). 

Finally, to evaluate the adaptability of the radiative 
package to a general circulation model, three heat bud- 
gets were attempted. The budgets were for the earth-at- 
mosphere system, the atmospheric column and the surface. 
Since heat transfer calculations due to evaporation and 
sensible heat are more easily made over ocean surfaces, 
London's atmospheric data was fitted to climatological 
ocean surface temperature and wind data from the United 
States Weather Bureau's Atlas of Climatic Charts of the 
Oceans (1938). Surface layer evaporation and sensible 
heat flux computations were not conducted poleward of 60N 
because of the possibility of ice cover. Rather, com- 
bined evaporation and sensible heat rates were taken from 


Vowinckel and Taylor (1964). 
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Figure l. Representation of the five-layer sounding 


model used for radiative transfer calculations. 


perature, mixing ratio, 


at the appropriate o-levels and k-levels (k 


Tem- 


pressure-scaled water-vapor 
(Au) and CO. (Ac) absorber mass increments are shown 
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100). 





II. DATA PREPARATION 


TEMPERATURE AND MOISTURE TERMS 

Northern hemisphere vertical and latitudinal distribu- 
tions of pressure, temperature, and relative humidity were 
taken from the climatological soundings of London (1957). 
The data was compiled in four mean seasonal tables. The 
latitudinal variation consisted of 10 degree latitude bands 
measured from the equator northward. In the vertical, data 
was distributed at 1 km intervals from the surface to the 
base of the stratosphere. 

With actual radiosonde data, e n information 
would be available above the tropopause to at least the 
0.1 sigma level. The G = 0.1 level was the highest level 
at which the water-vapor mixing ratio was considered to be 
recoverable from the fields of q. London's radiosonde data 
was continuous only to the tropopause. Since the tropo- 
pause varies in height latitudinally as well as seasonally, 
adjustment of the Ки distribution of T and q had to 
be made to insure that the data would be available to O - 
0.1 level. The simplest solution for the stratosphere was 
to assume isothermal conditions with constant relative 
humidity to the o = 0.1 level for any vertical sounding 
that ended below G = 0.1. 

It was necessary to interpolate from the distribution 


at successive km levels in the vertical (as listed in 


11 





London, 1957) to the specific sigma levels of Fipure 1. 
Temperatures were carried at the 1.0, 0.8, 0.6, 0.4, 0.2 
and 0.1 o-levels. Since surface data was available in 
London's soundings, Tig Was set equal to the listed sur- 
face temperature rather than requiring an extrapolation 


of the type 


Tig = 0.5 (3T, = T.) (1) 


Temperatures at other sigma surfaces, say at level k, were 
interpolated using a three-point Lagrangian formula of the 


form 


T (k) 


(T,) (P, Duell (Baba Bal 


+ 


(TD (PP 9) Ey PUR] -P,) (P,-P,) 


+ (TPL. (PP IP] PSI (PP) (2) 
where it is assumed that the data level (Po? To? lies above 
the o-level in question and (P,,T,) lies below that level. 

Mixing ratio values were calculated with London's data 


using 


4 = RH (0.622e )/P (3) 


The q-values were then interpolated to 0,-levels using the 


k 
Lagrangian formulation of (2) applied to q's at known 
levels. The q, values were centered at the odd O-levels 
of 0.9, 0.7, 0.5 and 0.3. These values are then considered 
to be the layer means between the even o-levels. 

Accurate humidity information is in general not avail- 
able above G = 0.3, especially in the stratosphere. To 
find a mean q-value representative of the layer (0,2) a 


different formulation was used. 


Je» 





Smith (1966) has shown a strong correlation between 


q(P) and the surface value 410 through the extrapolative 


formula 
442) (N (4) 
110 10 


London's soundings were tested to obtain a best-fit be- 


tween ln q and ln P of the form 


go (DA (5) 
43 = 
which is similar in form to that suggested by Smith. The 


HB ues of q(P) resulting from (5) for G = 0.9, 0.7, 0.5 
and 0.3 were used to solve for the profile parameter \. 
This parameter varied with latitude and season as shown 
in Table I. The relationship (5), with appropriate 
A-values, was then used to extrapolate q at 9 = 0,1. The 


resulting -values were then considered the mean 


= 
q-values for each layer (0,2). 
The procedure just described to derive q; was re- 
garded as justifiable in view of the following reasons: 
i) The residuals [q(est) - q(obs)], were small for 
11 k "and | 


(ii) The extrapolated corresponds to the proper 


. 
order of magnitude as listed in model atmospheres (Valley, 


1965) at pressures P = 100 Ы, 
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TABLE I. 


parameter À of Equation (5). 


Latitude 
0-10 
10-20 
20-30 
30-40 
40-50 
50-60 
60-70 
70-80 


80-90 


Winter 
3.1438 
SES 
3. 1205 
2.6437 
3.0002 
227210 
2555: 
Həz 


1.1743 


Spring 
2501768 
3.3024 
3.3978 
3.3838 
55176 
5::: 
2.7478 
25175 


176 


14 


Summer 
3.2885 
3.4381 
3.4272 
3.3563 
3.2480 
5551 
3.0684 
2.7886 


2.8434 


Seasonal and latitudinal values of the profile 


Fall 

7701202 
5. 2555 
3.2374 
3.0873 
3.1259 
3.1846 
5v 
2.4617 


2.1262 





B. CLOUD MODELING 

In both the solar and terrestrial radiative transfer 
processes, it is essential to know the fractional cloud 
Cover (CL) at each grid point. Clouds act as almost perfect 
blackbody radiators for infrared radiation, and as efficent 
reflectors for solar radiation. | 

There is a difficulty in parameterizing different 
cloud types and heights. In this model, as in the UCLA 
general circulation model, only a single cloud layer was 
parameterized for the purposes of inclusion in the vertical 
radiative transfer model at gridpoints. As is done in 
most other numerical prediction schemes involving cloud 
amount, CL, use is made of Smagorinsky's (1960) algorithm 
(applicable at approximately 700 mb) 

CL = 2,0 (een) - 0.7 (6) 

where 0 < CL < 1.0. When a cloud amount CL > O existed, 
it was assumed to extend between O = 0.8 to O =0,4. 

Comparison of the calculated values of CL computed 
from (6) and those tabulated from seasonal climatological 
estimates by Telegadas and London (1954) and Seinde (1954) 
showed a marked difference in the mid-latitudes as seen 
by comparison of Tables II and III. Since other climato- 
logical estimates of cloud cover, namely Sellers (1965), 
compared best with the data tabulated in Table IL,climato- 
logy was used in this model rather than the results of 


fable III. 
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TABLE TI. the seasonal distribution of total cloudiness 
(percent cloud cover) tabulated by London (1957). 


Lat N. Winter Spring Summer Fall 
0-10 47 əl 54 53 
10-20 36 42 49 48 
20—30 38 42 42 41 
30-40 50 52 41 46 
40-50 59 59 55 56 
50-60 63 62 63 66 
60-70 58 60 66 70 
70-80 47 29 69 70 
80-90 40 55 64 60 


l6 





£ "TE IIL.İThe seasonal distribution of total cloudiness 
(percent of sky cover) calculated by the Smagorinsky 
algorithm. 


Lat N Winter Spring Summer Fall 
0-10 30 40 56 46 
10-20 qo 20 . 48 36 
20-30 10 10 30 20 
30-40 18 14 22 16 
40-50 34 26 30 26 
50-60 42 34 42 42 
60-70 46 40 50 46 
70-80 46 38 50 48 
80-90 34 28 50 40 


17 





Even though London's cloud-cover climatology was 
accepted as valid here, it should be realized that London's 
results included overlapping clouds at multiple levels, 
which was impossible with the parameterization using (6). 
An improved cloud model permitting cloud existence at 
two or more levels in the vertical should be attempted 
in future radiative heating systems. At the present time, 
the added complexity of the radiative physics is considered 
too costly in computer time for inclusion in the present 
model. However, it should be noted that Arakawa et al. 
(1972) have already developed the concept of a three-layer 
radiative cloud model for a 1972 version of the UCLA 
general circulation model. 

Another reason for the choice of the climatological 
CL must be considered. With the use of climatology, 
transient synoptic scale perturbations and disturbances 
in the soundings are smoothed or EE Therefore, 
the day-to-day variations of relative humidity tend to Mes 
smoothed out during the season. The resulting smoothed 
Native humidity distributions thus tended to be less 
than their normal peak-values and give too small a magni- 


tude of CL when employing the Smagorinsky algorithm. 


C. -ABSORBER-MASS DEFINITIONS: PRESSURE SCALED 
Since mixing ratio values have now been computed at 
the five odd o-levels (Figure 1), it is possible to deter- 


Mine the pressure-Scaled water-vapor absorber masses, 
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Au,» associated with the five corresponding layers. This 
is done in the model in such a way that the additive 
scheme 

u(6,10) = u(8,10) + u(8,6) TO) 
holds. The layer water vapor absorber mass, u(2k, 2k-2), 
has been centered at odd levels (2k-1). The scaled 


absorber mass is calculated using 


q AP P 
5: m 0072 


0 


vhich, as has been noted, is the absorber mass in the 
layer (2k, 2k-2). AP is a constant in the analysis scheme 


equal to 0.2T. The mixing ratio value, is the 


Dre 
mean centered value in the layer (2k, 2k-2). The pressure 


0, 72 


ratio factor (P is the scaling factor for 


2k-17”o? 
collisonal line-broadening after Moller and Raschke (1964). 
This latter formulation disposes of the necessity of 
temperature scaling in the usual Lorentz line-broadening 


formula (c.f., Danard, 1969), vhere the pressure-ratio 


scaling is linear and a temperature-ratio scaling factor 


is indicated by E) times the N.T.P. mass within 
a layer. The subscrípt "O" denotes an N.T.P. condition 
where Po = 1013.25 mb, To = 273.16°K, and Po is the N.T.P. 


value of density. 





If (8) is integrated from Pok = (0.2k) Po to P = T, 
the result is 
P n=k>1 
K Po) “b s 552857... 
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where by definition 


G = Рӱр/т -——- 


k (10) 


For с0,, Moller and Raschke (1964) proposed pressure- 
scaling similar to that of water-vapor but with a pressure- 
ratio exponent of 0.65. For simplicity, a single 
pressure-scaling exponent of 0.72 was utilized for both 


cases. The form of the layer-reduced absorber mass of CO 


2 
a, T.P. nene in the layer (2k, 2k-2) is 
-4 Ro 2 
AC51..1 = 3.14 x 10 (12/59 ¿E CIT) 
The process of integration from Pok to Po yields the 
formulation 
ВЕ ор" 505“): RE 
“2k r? z 
where H = 7.995 x Oe en is the height of the homogeneous 


atmosphere in N.T.P. cms. Adjusting (12) to a O-coordinate 


system with P = T at the surface yields 


BIS x 107 


' _ 177 
Sok m” 


0 (13) 


Hİ em S ape 7 
0 


Table IVshows the seasonal distribution of the total 
pressure-scaled water-vapor absorber masses, u(0,10). As 
expected, the summer season contains the largest absorber 
mass with the least in winter. Spring and fall values are 
intermediate between the maximum and minimum values.  Com- 


putations for the CO, absorber-masses are dependent upon 


2 
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TABLE IV. Seasonal distributions of total pressure-scaled 
water-vapor masses, u(0,10), in gm cm?. 


Latitude Spring Summer Fall Winter 
0-10 777: 25777 4.023 34028 
10-20 2:2/1 4.316 ` 3.826 2.789 
0 2.417 3.601 35 055 hı 
30-40 dh ee) 2.897 2.133 77770 
40-50 1:099 5772 1.426 0.747 
50-60 00702 J. 795 0.932 0.449 
60-70 0.437 1.428 0.598 0.256 
70-80 05757 1.045 0.471 0.163 
50-90 0.192 0.901 0.319 0.105 
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the total mass of air in a unit column, and hence on the 
surface pressure. Since the variation in surface pressure 


averaged across all latitude bands was less than 10 mb, 


Нә 


the co, absorber-masses were nearly constant [c(0,10) 


125 cm/cm”). 
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LLL TELEREPTKLAL RADIATIVE TRANSFER 


A. THEORETICAL AND EMPIRICAL BACKGROUND 

Virtually all of the radiation packages in present use 
with large scale circulation models in the United States 
[cf., Manabe and Strickler (1964), Oliger et al. (1970), 
Arakawa et al. (1969)] make use of Yamamoto's (1952) water- 
vapor flux emissivities as a basis for their long-wave 
radiative flux calculations. The water-vapor flux emissi- 
vity is normally thought of as a function of both water- 
vapor scaled mass and temperature. Yet, in the normal range 
of temperatures (T > 220K), the emissivity is virtually 
independent of temperature. At temperatures T < 220K, the 


wave-averaged flux emissivity is defined after Yamamoto 


(1952) by 
co dB. 
f [1-1,, CL, 00] qr (17220K) dv 
pu T) = _— _ _ —_—_—_——  _—_ __ (14) 
/ CH 
^ aT (U,T) dv 


flere che numerator has its temperature dependent terms 
fixed at 220K,while the denominator equals ASt (T?) even 
when T<220K. The temperature affects the population of 
different rotational energy levels within the so-called 
"rotational" band. At low temperatures, Planck's function 
shifts maximal radiant energy to the longer wavelengths 


in the rotational band, according to Wien's displacement 
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law of 


À = 2897/T (15) 


max 

where 5. is in microns. Thus with T < 220K, the highly 
absorbed rotational lines are the only contributors to 
water-vapor absorptivity. This leads to an increase of 
absorption near the center of the rotational band and thus 
to an increase in mean absorptivity with decreasing tem- 
perature. This phenomenon is depicted graphically in 
Figure 2 after Sasamori (1968). 

Yamamoto's radiation chart has been the basis of most 
recent radiation calculations. Sasamori (1968, 1970) used 
Yamamoto's eae š theory to write empirical formulas 
for all of the emissivities needed for long-wave radiative 
flux calculations. These empirical formulas have been 
shown by Sasamori to yield adequate accuracy and were 
adapted for this numerical study. The temperature indepen- 
dent emissivity for water-vapor, from SER. is given by 

€(u) € 0.24 1081 (о + 0.010) + 0.622 (16) | 
25  tormulation vas used for T > 220K. For T < 220K, 

77 “nori gives the temperature-dependent part of the emis- 
sivity function by vater vapor as 


0.35310g4,g9u-0.44 -0.345510g,gu-.705 


e(u,T)=(8.34T Ju (17) 


When absorption by CO, and water-vapor are considered 


2 
jointly, the product law of transmissivities must be 


employed. In doing this, the transmissivity of the two 
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EMPIRICAL FORMULA 
« -40 C 
e -60 C 
e -80'C 


€ (u) 





Figure 2. Temperature dependent part of the emissivity 
with varied amounts of scaled water-vapor. 
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absorbers acting jointly may be expressed as 


тӱбмӱс) = трба)туСе) (18) 


The flux integral as formulated by Yamamoto is 


( 
Fe f Ba dB (19) 


Bees 


where u(T) is the water-vapor sounding mass measured away 
from the reference level. Equation (19) then becomes with 


the inclusion of co, 


u,c (T) 
P= [|  (li-rp(u)t,(e)ldB (20) 
u=0,6=0 
or 
u,c 
F = Ip [e€ (u,T) + Ag] dB (21) 
u=0,C=0 


The Ag term represents the effective emission due to co, 


not already counted in the water-vapor flux calculations. 
The added emission is expressed on the Yamamoto chart by 
an increase of emissivity at all points of the radiative 
sounding (except at the reference level points). The Ae 
term includes only the 15um band emitted flux of CO, vhich 
has not already been counted in emission by overlap vith 
the rotational band of vater vapor. Sasamori considered 
this co, band emissivity to be temperature independent 


in developing the empirical formulation 


10” 077 0 с)  ||1оре41.0641 (22) 
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In the following flux computations, the emissivity, €, 
will include both the effects of CO, 2176u7577(77) and that 


of water vapor given in (16) and (17). 


177 TECHNIQUES OF NET FLUX COMPUTATIONS 

Schematic representations of the Yamamoto radiation 
chart were utilized for considering the integration of the 
net long wave flux calculations. In Figure 3, the hatched 
area is a representation of downward flux at the surface. 
The area may be computed using the trapezoidal summation 
rule as follows: 

FY(10) = 0.5(16 (8-10) (8) 0-52) 1:16 (8-10)8€ (6-10)) 

(Bg-Bg)*[e(6-10)*c(4-10) ] (8,78, ) c [e (4-10) 


+ e(2-10)](B,-B,)+[e (1-10)+ e(2-10)] 


B(T 
(B,-B,)* 207 


B= 0 


ER 
e[(0-10),T]db) 


(23) 


Here B = St * T“ is the Stefan-Boltzmann blackbody flux 
arid e(8-10), for example, represents the emissivity com- 
puted from Sasamori's empirical equations using the water 
vapor and co, present in the layer (8-10). The term 

£ = €(u,T) + Ag iis solved by g(u,T)given by(17)and Ae by 
(22). The subscripts on the Stefan-Boltzmann fluxes (i.e., 
Bro” B 2) are identifers for the layer contribution from 


which the flux-emission takes place. 
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Figure 3. Schematic area-depiction of net flux at 
level 10. 





B(T,) D 
ii final integral of (23), 7 € dB , 
B=0 


may be evaluated by the integrated result 


ДИ ны 8.0 
4 = =— —— n r — — — 
Pt Er tela T1) (57335 log u + 3.56? S 
700 
0.26 
0.07262[1-0.62556(u + 0.286) J [log] ocot 150671) 


(24) 
In (24), it is understood that the integration spans 
the temperature range (0, T4) along the final isopleths 


Uo and ES @ e., Figure 3). The net flux at level #10 is 


given by 
* 
Fig = Bio - FY(10) | (25) 
where the asterisk superscript denotes F(up) - F(down) at 


the level. 
The net flux at level 6 is represented schematically 


by minus the hatched area in Figure 4. The numerical 


BIO 
modeling follows (23) where 


Bə = Big 7 0:5 (6(6-8)(Вӱ- В,) * (6(6-8)4 6(6-10)1 


(B107 Bo) + E (4-6) (B¿- By) + [e(4-6)+ e(2-6)] 


(B B,) + [e(2-6) + e (1-6) ](B,- B1) + 


.. 


B(T ) _ l 
NES 17076) Tİ dB) (26) 
B=0 | 
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BD = 


Figure á, Schematic area-depiction of net flux at 
level 6. ` 


/” 
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1000107 S rhe schematic for net flux at level 2. 


Following the previous computations 


Ju " Bio - 0.5 (62-4) (5,-B,) [e (2-4). €(2-6)] 


(Be - B,)+le(2-6)+ € (2-8) ] (B,- B.) + 


[e (2-8)+ SE EH 55 p 
B (T) _ 

+2 f € [(0-2), T] dB} 227) 
B-0 


* * 


* 
F F are clear sky 


The foregoing net fluxes Fio > Fo > Fo 


calculations only. A modification due to cloudy skies is 
made. The cloud-covered fraction is given by CL and the 
clear sky fraction is (1-CL). The cloud-covered fraction 
is considered an overcast which extends vertically from 
level 8 to ievel 4. Figure 6 shows a schematic flux dia- 
gram, after Yamamoto (1952), of net fluxes at levels 10 
and 2, when only the overcast sky fraction is considered. 
The horizontally hatched area is now the net flux at level 
10. The net flux at level 6 is zero since a layer of ap- 
proximately 50 meters in depth with average cloud droplet 
content may be TERIS. as a blackbody (Brunt, 1939) at 
the mean temperature of the cloud layer. The net flux at 
level 2 is the slant hatched area in the diagram. 

The net flux at level 10, with overcast skies is given 


by 
k 


Fio 7 Bio 


_ Bo) Q ” 0:56(6-10)) (28) 


Normally (28) represents a cooling effect at the earth's 


surface. At the same time, the net flux at level 8, or the 
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b 2 


B(T) —— 


Figure 5. Schematic area-depiction of net flux at 
level 2. 
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Figure 6. Schematic area-depiction of net flux at 
reference levels 10 and 2 when an overcast cloud 
extends between levels 8 and 4. 





cloud base, is an equal magnitude (Fig. 6) so that the 
layer (8-10) has zero flux divergence. However, the cloud 


base receives a net flux warming equal in magnitude to Fig 


The net flux at level 2 for the overcast case is given 


by 


F - B, - 0.516(2-4)(B,- B,)+ EE B.) 


2 4 1 


2:0 
+2 f 217022) 7 dE 
: . B=0 (29) 


* 
TaS net flux, F, ; is much smaller than that for clear 


skies because of the lead term B, on the right side of 


27577 as contrasted with B in (27): 


10 
The net fluxes, computed at the three reference levels, 


must now be cloud-weighted to yield final results for net 


long-wave flux transfer across the interface by 


* * de 
Fig = 17/00 (7521 CL 4 [E (252 1(1-CL) (30) 
рд" = [FÍ(26)] (1 - CL) : (31) 
Fo. = [F5(29)] cL + [*(27)]@-cL) 7995 


The temperature-dependent flux integral denoted by 
c (u: Ti) dB was evaluated for an entire range 
of final Tj and u, values which 


could occur in the atmosphere. These results are listed 


in Appendix Table A-l. 


er DIVERGENCE OF NET FLUX AND RESULTANT COOLING RATES 
The resultant, cloud weighted, seasonal distributions 


of net flux at the three reference levels of o = 0.2, 0.6, 
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and 1.0 are listed in Tables V(a-d). The 24-hour tempera- 
ture changes between the reference levels due to terrestrial 
radiation alone are also listed. The mean radiative tem- 
perature change centered at the midpoint of the column was 
computed by the first law of thermodynamics so that the cool- 


ing rate for the layer (6,10) is 


oT 1000 








EE EI 3 
25277) 01” e, (33) 
and for the layer (2,6) is 
ƏT _ * 1000 
= (E, F. 7 070 5 ES > (34) 


where T is surface pressure in mb. 

To compare the accuracy of the long-wave calculations, 
the meridionally averaged mean values for the fluxes at the 
levels o = 0.2 and 1.0 were computed. The result was repre- 
sented by the percentage of the meridionally averaged inso- 
lation at the top of the troposphere. Katayama's 63 percent 
(1967), and London's 62 percent (1957) for net flux at 200mb 
are in good agreement with the 64 percent from this study. 
The net flux at the surface was 18 percent for all three 
studies. Converting these spatially averaged results to 
annual means (ly da ve get a net outgoing radiation at 
P * 200 mb of 430 from the Katayama-London means compared 
to oux 450. At the surface, the net outgoing long-wave radia- 
Bron is 120 (ly 152 compared to our 131. 

In the net flux calculations, a comparison vas made be- 
tveen three different atmospheric um ons in order to 
evaluate the contribution of cloudiness to the total outgoing 


flux from the earth-atmosphere system. One test considered 
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only clear skies, another had average cloudiness and the 
third tested totally the overcast model. Results were com- 
piled in mean annual latitudinal distributions of long- 
wave net flux in Tables VI(a-c). 

Since temperatures normally decrease with height in 
the atmosphere, a positive net flux occurred at all levels. 
Strong Surface temperature inversions can lead to negative 
gə but these conditions were not present in this study. 
Also, in an atmosphere where the effective optical depth 
decreased sharply with height, the net radiative flux in- 
creases sharply with height at the levels of the decreased 
water-vapor mass. This happens primarily because most of 
the emitted flux from the lower layer passes upwards  un- 
attenuated by the dry upper layers. 

The net flux through the tropopause is larger at all 
latitudes for clear skies than for conditions of average 
or overcast cloudiness. This occurs because the net outgo- 
ing radiation comes from a lower, and therefore warmer 
layer in the case of clear skies. In the case of overcast 
skies, the radiating layer is the cloud top at Oo = 0.4 
Rather than о = 1.0 for clear skies. The tropopause net 
flux with average cloudiness is a weighted average between 
the overcast and clear-sky cases of Table VI. 

The largest variability of net flux magnitudes 2301 cloud 
cover occurred at the G = 0.6 level. The minimum value of 
zero is observed during total cloudiness since this level lies 
within the cloud layer which was treated as a blackbody. Maxi- 


mum values occurred of course in conditions of clear skies. 
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IVa SOLAR RADIATION TREATMENT 


AM EFFECTIVE SOLAR INSOLATION 

The solar radiation treatment is not as straightforward. 
Following Joseph (1966), a solar constant of 2.00 ly SS 
was assumed. A 4% attenuation by oxygen and ozone above 
the tropopause was assumed which gave an effective solar 
constant of 1.92 ly 250 The effective tropopause height 
was assumed to be at G = 0.2. 


The effective insolation at the tropopause was computed 


from 
F = s(r/r ) ° D cos Ü (35) 

where 

S = effective solar constant 

r/r = radius vector of earth relative to the sun 

8 = zenith angle of sun relative to observing point 

cos 0- daily average of cos 909 

D = fractional length of daytime 


The (r/r |) is listed as a function of the Julian date 
in the ieh sonian Meteorological Tables (List, 1958). The 
following table, after Manabe and Moller (1961) was employed 
to derive the daily mean ək O and duration of insolation 
as a function of latitude and time of the year.  Computa- 
tions of F were based upon mean mid-seasonal values of 
Spring ~ April 15th 


Summer - July 15th 
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TABLE VII. 


Lat. 


15 
2. 
35 
4 5 
55 
65 
22 


85 


Fractional Length 


of Daytime 


Apr 
.508 
5-1 
2513 J 
546 
2502 
. 596 
629 
7750 


1.000 


Jul 


.517 


5 


562 


.596 


.637 


. 708 


55 


1.000 


1.000 


Oct 


«200 
492 
WER 
2.71 
454 
A37 
.404 


2329 


Jan 
. 496 
~471 
.450 
.421 
762 
5: 


. 208 


43 


Weighted Mean 
cos 0 


Apr 
2025 
. 618 
5: 
5556 
2201 
.423 
e AN 
.241 


.168 


Jul 
.587 
.601 
57 
.567 
2523 
453 
. 369 
. 311 


. 318 


Oct 


.614 


ASIE 


.524 


.458 


. 379 


25 


.176 


. 071 


ln. con 17 ond latitudinal distributions of 
duration of insolation and daily mean values of cos Ö, 


Jan 
‚591 
(549 
(474 


393 


-.17 


“05 


.106 





Monee Vill. Incoming solar radiation on a horizontal 
surface at the top of the tropopause (ly day71), 


Latitude Spring Summer Fall Winter 
0-10 884 867 844 784 
10-20 896 5 783 692 
20-30 889 952 696 571 
30-40 848 965 5 443 
40-50 784 948 473 306 
50-60 702 916 7: 174 
60-70 604 882 195 59 
70-80 503 888 6.4 0 
80-90 468 908 0 0 


4 4 





w 


PPP 





"9/1/90 








Fall - October 15th 
Winter - January 15th 
Joseph's method (1966) was followed in partitioning the 
effective solar insolation into a part subject to water- 
vapor absorption but not to Rayleigh scattering, and the 
remainder subject to only Rayleigh era The presence 
of clouds introduced cloud reflectivities im both portions. 
Thus, the solar insolation subject to tropospheric absorp- 
tion comprises wavelengths > 0.9 um and was obtained from 
..: (36) 
In the vavelengths < 0.9 um, water vapor absorption is 
negligible and the solar insolation subject to Rayleigh 


Scattering was obtained from 


Fo = 0.651F (37) 


DISPOSITION OF FLUX SUBJECT TO SCATTERING 

In treating the fraction of the solar insolation subject 
to scattering, Joseph (1966), found that the Rayleigh 
scattering values at sea level in clear skies (after Coulson, 
1959) could be fitted by least squares to the form 


| 38 
ас = 0.085 + 0.25074 log E ess 0) DER 
o 


where T is the observed surface pressure. For a full-day 
the value of sec 8 in (38) was treated as sec 9, the time- 
mean sec 8, which was taken as the inverse of the mean 


eos O values from Table VII. 
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Another reflective parameter is the surface albedo, a 
After numerical experimentation designed to test for a 
global radiative heat balance, a constant average surface 
albedo of — - 0.14 vas deduced. It should be noted that 
this value of = turned out to be equal to that deduced by 
500050 (1955). A discussion of the validity of a constant 
surface albedo will be reserved for the surface balance 
section. 

Considering the possibility of multiple reflections be- 
tween earth and atmosphere, each of which diverts downward 
the fraction o of the earth surface reflectance : the 
clear sky insolation at earth after scattering is given by 


Iso 


Ы o a t--- (aga) + .5: (39) 


Š 
or 


Iso 


Fs (1-6) (1-a J/ (1-00) (40) 


For the existence of a cloud cover, the Rayleigh  re- 
mieetion coefficient of (38) was applied at the cloud top 
zi 0.0, rather than 0 = 1. At the cloud top, a unique 
cloud reflectivity fraction R of 0.5 was chosen everywhere. 
This reflection is basically a Mie-type, diffuse back- 
Scattering caused by water droplets in the clouds. 


Thus, at the cloud top, the entering solar-beam insola- 





tion is given by NT 
Fo [l-a po sec 9)](1-R_) 


Je en s l-a R (41) 
5 C 


This downward scattered beam is forced into a diffuse 


angular array by the Mie scattering process. In order to 
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25 “int for Turther attenuation by Rayleigh scattering vhen 
the ray paths are diffuse, a mean slant path angle 9 such 
Haut sec 95 = 5/3 was judged more appropriate than the actual 
755: (Katayama, 1966). 

The Rayleigh scattering coefficient beneath the cloud 
top (i.e. in layer (4,10)) was found to be 


72 “sp, = “sp, - 0.25074 log 2.5 (42) 


Thus, be lowernicectomcetop has a constant value of 0.0998. 


Zen 


In this connection, it may be noted that Manabe and Strick- 
ler (1964) employed a constant value of a. = 0.07 S cos 8 
loss to space by Rayleigh scattering regardless of the pres- 
ence of clouds. 

Multiple reflections take place below the cloud top. 
Martin (1972) has shown that upon simplification, the inso- 
lation at earth after cloudy-sky scattering is 


Fsca (7sp77y (177sp? 10-da) (1-a IR, 


lec" BU ten (1-8. ) + 


S 


(1-Е, (4-0, 2 


(43) 
In the above accounting, certain small atmospheric re- 
flection terms occurring in the direction cloud to earth 
(or reverse) have been ignored as small compared to the air 
transmission factor (1 - Ron)‘ 
Finally, the weighted mean earth-surface absorption of 
the 20: solar beam becomes 
LI. (10) = (1-CL)Igg SR Ig (CL) (44) 


where I is given in (40) and I in (43): 


S 0 S C 
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TROPOSPHERIC ABSORPTION OF SOLAR RADIATION 


In considering the absorption of the F, component of 


A 


solar insolation, the Manabe-Moller solar-absorptivity 


function 


0.303 


a = 0.0946 (u sec 0) (45) 


22K 
was used. The form of Equation 45 indicates that the ab- 
Sorption function a is applied over the pressure-scaled 
water vapor mass in the layer (2,2K) along the zenith 
: path of slant angle 0. Restricting absorption to the ratio 
0.349F requires that the multiplicative constant 0.0946 of 
155) be adjusted to (0.0946/0.349) = 0.271. 

The absorbed insolation in the layer (2,6) is then 
0. 303 


A(2,6) = 0.271[u(2,6)sec 6] Fa (46) 


Similarly, the absorbed insolation in the layer (2,10) is 


0.303, 


A(2,10)= 0.271[u(2,10)sec 9] E 


(47) 


Therefore, the absorbed insolation in the layer (6,10) is 
given by 

771075 5:60 10) - A(2.6) (48) 
Since the amount of solar beam absorbed in the layer (2,10) 
is available BV (АТ che solar insolation in the F,-parti- 


A 


emn received at earth is given by F,- A(2,10). This por- 


A 


tion is reduced by ground albedo E so that the 


earth-absorbed insolation resulting from FA is 


ə 0 
0.3 3) 


1 (10) = F,{1-0.271[u(2,10)sec 0] (1-0) (49) 
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modification is made for cloudy skies. The absorption 
in the layer (2,4) was computed according to the Manabe- 
Moller absorptivity function [Equation (46)]. The solar 


beam at level 0.4 is given by 


| F ACA) - FA Er) (50) 


Of the F “energy available at level O = 0.4, the fraction 
(1-Е,) is transmitted through the cloud top. Thus 


F, (4) = [F -A(2,4)](1-R_) (51) 


Within the cloud layer, the water vapor path below G = 0.4 
is subject to the mean diffuse-path augmentation factor 5/3 
(after Katayama, 1966). 

Inside the cloud layer, another variable has been intro- 
duced. The value is u or equivalent water vapor mass for 
droplet absorption in clouds. The values of the absorption 
coefficient aCLD due solely to cloud particles have been 
taken after theoretical computations of Korb, Michalovski 


and Moller (1956). These values are listed after Manabe 


and Strickler (1964, p. 367) in Table IX, below. 


TABLE IX. Absorptivities (aCLD” of cloud droplets. 


Latitude Band as LD 
0-39N, low-to-middle clouds 0.035 
40-64N, low-to-middle clouds 0.020 
65-90N, ice-crystal clouds 0.005 


Since only a single cloud layer was used in the radia- 
tive model, the larger absorptivities associated with the 


larger water content was identified with the 0-39N latitude 
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band of Table IX. The values of acıb are due solely to ab- 


sorption by either liquid water or the solid ice content of 


clouds. The values of as LD are considered mean values for 


the latitude bands shown in Table IX. Moreover, the LD 


values of Table IX are averaged for the entire solar beam, 


and not merely for the F, portion of the insolation. Hence, 


A 


the value of 8с11) has been multiplied by (1.0/0.349) for 


£5— ocoation to the F -partition. In this spectral region, 


A 


the effective cloud-vater absorptivities become 


* 
a = 2.8653 2 CLD (52) 


or 


* 


a^- (2(u(4,8)*u,) ] Ep 


0.303. 


[2u(4,8)] 171 (53) 


| * 
Equation (53) may be solved for u as a function of a and 


u(4,8) 
R : 0067790 01958 
с a ° ° ә 
и(4,8) Cilt OP ca 8) J -1.0) 29) 


Values of u. were computed as a function of u(4,8) for the. 
latitudes of Table IX and are shown in Appendix Table B. 


With u. known, the following calculations are made 


1 u 
17 6)- Р. (4)(0.27112 CA oJ IÓ (55) 
A(6,10)2F, (4) (0.271[2(u (4, 10) *u ) 19 * 772) -A (4,6) cə 


vhere F ACÓ) is evaluated upon entering the cloud top. 
Finally, the insolation reaching earth after absorption in 


cloudy skies is 


Tic (10) ED - A(4,6) - A(6,10) (57) 
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10270 the Surface albedo effect 55 included, the earth- 


absorbed part of O) becomes 


Poet O Eesen) (58) 


D. NET INSOLATION ABSORBED AT EARTH 
The weighted sum for the insolation at earth after ab- 
sorption is given by 


Т.с̧(10) = CL[I c 0.55 EE (1-CL)[I,g ae (59) 


Thus the total insolation absorbed at earth is given by 


205 1: 10077 (10) (60) 


The resulting values of Q10 are listed in Table X. 
E. COMPARISON OF MERIDIONALLY AVERAGED RADIATION 
QUANTITIES 
In order to compare the various annual outputs of this 
radiation model with the results of other researchers, the 
radiational quantities derived by the NPS model were time- 
averaged and space-averaged over the Northern Hemisphere. 
The resulting global mean annual quantities are compared 
with the results by Katayama (1967), Budyko (1956) and 
London (1957). The mean values are presented as a percent- 
age of the solar radiation at the top of the troposphere. 
Meridionally-averaged values were generated using the 


formulation 


s m... (61) 
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TELE X, Seasonal distribution of insolation absorbed at 
earth LH: ol in 17 day-l 


Latitude Spring Summer Fall Winter 
0-10 451 430 424 409 
10-20 484 471 404 385 
20-30 437 >30) 374 316 
30-40 450 525 315 234 
40-50 405 484 241 157 
50-60 361 451 166 88 
60-70 517 430 95 30 
70-80 264 430 31 0 
80-90 248 454 0 0 
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where 


R = annual mean radiative quantity at latitude ф 

a = radius of earth 

b = center of the latitude band under consideration 
do = latitude increment in degrees between successive 


bands. This increment is a constant of 10° for 
this study. 


The formulation in Equation (61) may be reduced to a 


series summation using cosine weighting factors as indicated 


in 
9 
ñ 2 R cos $a 
E + (62) 
9 
> соз ф 
n=l 5 


where $1 = 5°, P = 15°, etc. 

When the results of this study are compared with those 
of Katayama, Budyko, and London in Table XI, it is found 
that they are in rather good agreement a a whole. However, 
ae are some minor differences. First, the older model 
estimates of solar radiation absorbed in the atmosphere are 
generally 3 to 5 units higher than determined by this study. 
Also, while the average global albedo, A (Ф), determined in 
this study is typical of usual values quoted, namely 34%, 
the meridional distribution of a was more nearly uniform 
with respect to Q than would have resulted if the means had 
been based upon hourly computations for the months involved. 
This result leads to a lower than normal value of solar 


radiation absorbed at the surface in low latitudes and a 
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TABLE XI. Comparison of meridionally averaged values of 
radiation quantities. The values are presented as per- 
centages of the global mean insolation at the top of the 
troposphere. 


Katayama Budyko London Plante 


Ear radiation at 


top of troposphere 100 100 100 100 
Solar radiation absorbed 

in atmosphere 18 16 14 13 
Global albedo 37 40 35 34 


Solar radiation absorbed 
at surface 45 44 47 52 


Outgoing long-wave 
radiation at 200 mb 63 60 62 64 


Outgoing long-wave 
radiation at surface 18 17 18 18 


Sensible and latent 
heat flux (see Section IX) 29 27 29 31 
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cer than typical value at the surface in high latitudes. 
memes reason for this is basically threefold: First, atmos- 
pheric absorption by dust was not considered in the NPS 
model, since it is essentially a predictive model. However, 
London (1957) applied the results of Houghton (1954) who 


used atmospheric transmissivity due to dust as given by 


m 


un (0.95) (63) 
where m = sec 0, represents the air mass traversed by the 
solar beam. London's depletion by atmospheric dust would 
then be given by 

D = (1 - T.) (64) 


of which 254 is considered subject to absorption. Even 
50, the contribution of dust absorption is small and can 
generally be considered of the order of lZ, which is 
approximately the difference between the NPS absorption 
figure and that of London's. 

Another contribution to absorption not previously con- 
E red here is secondary absorption of the reflected por- 
tion of the solar beam. Since the reflected portion of the 
solar beam above clouds is subject to the same percentage | 
absorption as the incoming beam, the error in absorption 
may be as high as 1% of the incoming solar radiation. Neither 
London's study nor this study considered secondary absorp- 
tion of the reflected beam and consequently the resulting 
atmospheric absorptions are closely comparable. Katayama 


considered secondary absorption following a first reflection 
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urbe solar beam by a cloud or by ground or both. Asa 
result Katayama's solar absorption in air was the largest 
of the four cases considered. 

ihe third consideration is Katayama's cloud cover. In 
the Katayama radiative study, a larger climatological cloud 
cover was considered than in the present study, or that of 
London or Budyko. The cloud MET were generally charac- 
teristic of coverages over the oceans (after McDonald, 1938) 
rather than a combined ocean-continent regime. 

It is not exactly clear how Katayama allowed McDonald's 
overlapping cloud-layers to be combined, although he con- 
strained the total cloud-coverage to agree with that of SIr 
Napier Shaw's (1936). The significant feature of the 
greater cloud coverages is that they afford not only greater 
solar absorptivity but also a greater global albedo as 
postulated by Katayama's empirical formulas 


A 


S Az [1-exp(-b, h) ] (65) 


R 


... R. [1-exp(-b, h) J (66) 


b A R are 


where h is the cloud thickness and bi 29 åo? Ro 


characteristic values for a given cloud type. Thus Kata- 
yama finds a greater cloud absorptivity than was found by 
London or Plante (Table XI) as well as a larger global 


albedo because of the larger cloud coverages. 


F. GLOBAL ALBEDO CONSIDERATIONS 
In the process of deriving the transmission of the two 


solar components F, and F the following dispositional 


A S 


terms were discussed: 
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(i) ane combined insolation rate at the tropopause 
@emoted by F of Equation (35). 

(ii) the absorption rate of insolation in the atmos- 
phere 

1772717 the absorption rate of insolation, Q10> at earth 

With these terms, the reflected E rate from 
planet earth, as a function of latitude and season may be 
computed. This rate may be denoted by 75 where SE is the 
global albedo. The term SEN must be identified as the 
difference of 


Daach A F (67) 


F - A(2,10) 

It should be noted that three significant reflective 
parameters have affected the value of the left side of (67) 
namely the sky albedo, the cloud albedo (R. = 0.5), and the 
earth-surface albedo, 5 For simplicity, the latter tvo 
vere held constant at reasonably vell accepted figures in- 
dependent of latitude. While direct computations of the 
right side of (67) were made at each step of computing (ii) 
and (iii), the resulting global albedo values must also be 
consistent with the e Uds inferred from the left side of 
567). 

A note of caution in regard to comparing global albedo 
as shown in Table XI. A correct mean global albedo may 
distort the atmospheric heat balance if not distributed 
correctly Een respect to latitude. For example, if a 


meridional bias exists in the computed global albedo, Ab 





with too large a value at low latitudes and too small a 
value at high latitudes, the result will influence the pole- 
ward heat transport required for balance, after independent 
Euputation of latent and sensible heat transfers across 

the earth-air interface (see Section IX). This type of 

bias in global albedo appears to have developed in this 
study presumably by oversimplifying the ray-paths in clouds. 
Choosing a single mean path through clouds for a given re- 
presentative day, and time of year seems to have overesti- 
mated the global albedo in this study in lower latitudes 


while underestimating it at high latitudes. 
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V. RADIATIVE TEMPERATURE CHANGES 


Long-wave radiation, with a normal temperature lapse, 
will lead to a cooling of the atmosphere as seen as seasonal 
distributions in Tables y(a-d). In certain locations and 
seasons, some of the atmospheric layers considered were 
subject to temperature inversions which led to local heating 
rather than the usual cooling. Another heating situation 
examined was at an overcast cloud base. The terminology, 
"cloud base warming' or 'green house effect' may be used 
to describe this phenomenon. Since the net flux within a 
cloud layer is zero, a negative net RU increment existed 
in the cloud to sub-cloud layer (8,10). With the sign 
convention adopted here, this represents a base-warming 
effect. In the cloud top region, there is a substantial 
mete tlux increment FQ = Fo which leads to considerable 
cloud-top cooling. Thus, in infrared radiative temperature 
changes the role of clouds ís clearly to increase the cooling 
rate at the cloud tops and to decrease the cooling rate 
at the cloud bases. This process can be considered as a 
decrease of the atmospheric stability within the cloud layer. 

The mean annual infrared temperature changes in the 
troposphere are plotted in Figure 7 as a function of latitude. 
The dependency of these cooling rates upon the climato- 


logical cloud cover at these same latitudes is evident by 
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inspecting Figure 8. A linear regression was performed re- 
lating the mean long-wave temperature changes of each tropo- 
sphere layer against the total cloud cover used in the 
computation at latitude db. Each latitude band for each 
season was considered separately, giving a total of 36 cases. 
For the layer (6,10) the correlation coefficient was -0.95 
while for the layer (2,6) the correlation was 40.99. Graphs 
showing the comparison of temperature changes found with 
clear and with average cloudy skies for the layer (2,6) and 
layer (6,10) are shown in Figures 9 and 10, and clearly 
demonstrate the phenomena of cloud-top cooling and cloud- 
base warming. I 

In conditions of normal cloudiness, the infrared heat 
losses are partially offset by atmospheric absorption of 


solar radiation. The resultant temperature change in the 


layer (6,10) is 
* * 
. - Fig ) + A(6,10) 


oT MER Co 109 NS E 
ae (5,10) a 0.47(1000/g) e, ` 


(68) 
The temperature change in the layer (2,6) is given by 
* * 
- (F, -= ur 030) 
0.47(1000/g) = 


oT 


тели (69) 


In most cases, this absorption of solar insolation is smal- 
ler than the cooling that results from long wave radiation 
and acts only to modify the pattern of radiative heat loss, 
as shown in Figure 11. The dashed curves represent the 
annual mean rate of temperature change due to long wave 


radiative cooling for the layers (6,10) and (2,6). The 


61 





о 


06 


*әрајтзет YITA UOTINQTAISTP 19A09 рпото [enuue ueoy 


30011171 
08 oz 09 os Or oç oz 


*д әлп8тд 


ol 





007:::o0 pnoT2 TE3OL 


62 





06 


*SITAS Ápno[o e2ge?aoA? pue i?Po[o 


u3TA (9*2) zəKer UT әјел 3SurToooS5 pəzezjur [enuue uegoy 


3Q n111V1 
08 OL 09 os ot 06 02 





..—— sonunun. 


Ku < fo” — 


ssəurpnoT2 ə3erəAV 


“6 әлпатд 


Ot 


(,, AVG Do) 39NVH? 33n1V334V31 


63 





06 


"soTys ApnoToO eSeioAe pue gar? 


437” (0159) AB9ÁABT UT eaeaz 3uTTOOƏO peaieaigur [enuue ueosy 


309N111V1 
02 OL 09 os ОР Ot 02 


"OT sandra 


ot 





(,- AV Do) 39NVHO 3ufhiviudidw3il 


64 





*(9*2) Pue (OT°9) Sise/ier 103 UOT3ETPELN 3əu pue 3uTTooəo рәделјир 10] 
esgueys sanze1aduaz oriouydsom3e Tenuue uesu jo uosf1eduo) "TT 575: 





3anill1V1 : 
06 08 OL 09 OS ot 06 02 01 
— Yk 
— n ara 
» * 
5” EO 
o o =f 
(019) xəZeq x 
- 
m 
zı 
> 
à 
с 
2) 
m 
15504 
On de : Q 
ke 
-— ə 0 s > 
SOL Q 
кк 
5. eee ac = m 
ni e^ 
€ 
(952) 194e7 25 O 
: — 
“o Е-о 


(, AV 


65 





sn Curves represent the result of net radiative cooling 
(combined effects of long and short wave radiation) in the 
Same layers. Atmospheric absorption clearly only modified 
the magnitude of radiative cooling while the shape of the 
curves was unchanged. By comparison of Figures 8 and 11, 
Men Correlations between net radiatie temperature changes 


5200 total cloud cover may once again be noted, 
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VI. RADIATIVE BALANCE AT THE TROPOPAUSE 


A radiative budget was calculated at the tropopause and 
cOmpared to the results of Houghton (1954). Houghton made 
a comparison of the net incoming solar radiation at the 
tropopause with the outgoing long-wave flux. He determined 
that insofar as the earth-atmosphere combined system is con- 
cerned, a column at about 35°N latitude experiences a mean 
radiative energy balance. Corresponding latitude belts of 
the system located equatorward of latitude 35N have a mean 
annual radiative energy surplus. Those latitude bands pole- 
ward of 35°N latitude experience a mean annual radiative 
energy deficit. These mean annual radiative surpluses and 
deficits impose requirements for horizontal heat transfer 
rates across latitude walls in both air and ocean. 

The net insolation at the tropopause (4) 15 given by 


EE DECHE (70) 


o 
vhere Q. is the incident insolation at the tropopause after 
a 5” on due to attenuation in the stratosphere, Qe 
is the reflected insolation returned to space. The ratio 
of Q / [1.04 Q,] is the measure of the earth atmosphere 
planetary albedo. Qa is given by 

200500 0 (71) 
vhere QA is the short wave radiation absorbed in the atmos- 


phere and Q10 is the total short wave radiation absorbed 


at earth, 
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Bie resultant Qn is a function of the Rayleigh back 
m ering, cloud reflectivity and surface albedo. The 
factors determíning Rayleigh scattering were discussed 
in Section IV. 

Treatment of cloud reflectivity was limited here by 
the use of a single cloud layer. A global average of 
about one-fourth of the incident 5717 radiation is re- 
flected to space by clouds. In this model a single cloud 
albedo 5 0.5 vas used. Sellers (1965) employed a range 
Of albedos for differing cloud types. A similar treatment 
using various cloud types with various albedos was adopted 
by London. On a mean annual basis, 0.5 appears to be a 
good approximation for the hemispheric cloud reflectivity, 
but on an individual gridpoint basis the cloud albedos 
adopted in this study may be an oversimplification. Future 
testing in this area is envisioned for the NPS predictive 
model. 

A single constant surface albedo of 0.14 was chosen 
after 577 incntation with the tropospheric energy balance. 
Budyko (1955) estimated the mean earth surface albedo as 
Mea while Katayama (1967) found a value of 0.13. But the 
albedo varies with the type of surface. The short-wave 
albedo is highest over fresh snow ( 90%) and is lowest over 
water and forests (approximately 5%). Since there are 
seasonal as well as annual variations in the extent of 


vegetation and snow cover, an attempt at strict evaluation 
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of ground albedo for these variations is at best difficult. 
Thus, for simplicity, a single mean ground albedo was em- 
ployed. 

The net outgoing long-wave radiation at the tropopause 
was given by EI since the 0.2 sigma level was adopted as 
Bu mean height of the tropopause. 

Eheretore, the net radiative flux across the tropopause 
is given by 

Q. = Q. ` F, (72) 

Results for seasonal variations in the net flux across the 
tropopause are listed in Tables XII(a-d). The mean annual 
2 ve input by latitude bands is listed in Table XIl1(e) 
along vith the cosine-latitude veighting factors. The re- 
sulting mean Northern Hemisphere radiative budget is zero 
as anticipated, since it is assumed that the hemisphere is 
in annual radiative balance. 

During the summer, the tropopause has a net radiative 
surplus throughout all latitude bands. As seen in Table 
XII(b), the incident insolation is largest during this sea- 
son and far exceeds the small increase in net long wave 
radiation due to warmer summer temperatures. The opposite 
case is true for winter. Low insolation rates produced a 
negative radiative budget at all latitude bands except the 
0-10° band. This particular latitude band has a positive 
radiative net surplus during the entire year. Spring and 
fall show transitions between the extremes of summer and 


vinter. 
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Effects of variations in cloud cover on the computed 
net radiative input at the tropopause can be seen inm Table 
XIII which lists spatial mean averages of the radiative 
quantities across the tropopause. It is seen that clouds 


exert an overall cooling effect in this model. 


TABLE XIII. Spatial mean annual net radiative input at 
the tropopause (Northern Hemisphere). 


CL = 0 CL = 1.0 Average CL 


Q. 694 694 694 
Q. 130 360 244 
0, 564 334 450 
p 528 373 450 
Net Tropopause 
Radiative Input 36 -39 0 
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VII. GLOBAL ALBEDO 


The earth's global albedo results from the upward re- 
flection by the ground, by clouds and by scattering by 
atmospheric molecules of solar radiation. Variations in 
computed planetary albedo are strongly dependent on total 
cloud cover. The dístribution with latitude and season 
of global albedo is listed in Table XIV. There is a pro- 
nounced minimum in the subtropics, which is a region of 
minimum cloudiness. The minimum shifts equatorward during 
the winter and poleward during the summer, following move- 
ment of the ITCZ. The increased albedo in high latitudes 
during all seasons results from an increase in cloudiness 
5571 latitude up to about 60 to 70N. 

Total global albedo when averaged over all latitudes 
and seasons is about 34 percent. The total global albedo 
shows little variation with season and agrees closely with 
London's results. Figure 12 shows a comparison between 
mean annual planetary albedo with latitude computed by 
London and the results of this study. Even though the 
spatially averaged annual global albedo of 34% agrees close- 
ly with London, there is disagreement in its distribution 
with latitude, a fact referred to in Section IV. 

In tropical latitudes, London's global albedo was less 
than that computed in this study, while in the northern 


latitudes, his values exceeded those computed here. Since 
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TABLE XIV. Ses sonaldistribution of planetary albedo with 
latitude (in percent of 0). 


Latitude Spring Summer Fall Winter 
0-10 35 37 36 34 
1020 32 35 34 30 
Gan 32 32 SE 31 
30-40 36 32 34 36 
40-50 39 37 38 39 
50-60 40 40 41 41 
60-70 39 40 42 40 
70-80 40 41 41 5 
80-90 40 40 : = 
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muemtcical total cloud cover amounts were used in both 
studies, the treatment of atmospheric and cloud albedos 
seem to be the major cause for this difference. The 
present study used a constant cloud albedo of R. = 0.5 
while London used varying cloud albedos associated with 
constituent cloud types after Haurwitz (1948). London 
thus used a multi-layer cloud model while this was not 
done here. Furthermore, the device of replacing the 
diffuse path length in clouds by an equivalent path 5/3 
longer seems to give too short a path subject to cloud- 


reflecting processes at high latitudes and too long a 


ph at low latitudes. This is important if cloud-reflec- 


pance is to be treated by Mie-scattering theory. Varia- 


Eus of surface reflectance with solar zenith angle also 


plays a secondary role in explaining the variations of 


= with latitude. 
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Vitis ATMOSPHERIC T BUDCET 


A simplified annual atmospheric heat budget for the 
Northern Hemisphere, along with seasonal distributions of 
radiative surpluses and deficits was constructed. The 
Stratosphere was assumed to be in radiative equilibrium. 


The atmosphere radiation balance is shown schematically 


below: 
* 
Q. On E, 
E 1. a . мы 
X 
Qio d Q. E 


The net cooling contribution due to long wave radiation is 


given by 
* 
N O Си 


The effective insolation at the top of the troposphere 
is the total incident solar radiation minus the reflected 
insolation to space. The total atmospheric absorption is 
27 rore | 

Q, = (Q; - Q.) - Q4 (74) 
where Q10 is the short wave radiation absorbed at earth 
using a mean surface albedo, 2 = „14. A final heating 
term consists of the latent heat of evaporation and sen- 
sible heat transfer rates across the earth-air interface 
and denoted by Ор: 


50.4 (75) 
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50 Pparameterization of Qr and Qe Homcetat ted sin Section Ix. 
For each latitude band, the net atmospheric heating rate by 
the three mechanisms (insolation absorption, long-wave 
absorption, interface convecticn) is given by: 
5771200. (76) 
The resulting atmospheric heating (cooling) rates in 1y dc m 
are summarized in Tables XV(a-d) for each latitude band and 
each of the four seasons, while Table XV(e) gives the annual 
mean, derived by averaging the four seasonal heating rates, 
as a function of latitude. 

NIIS to be noted that Qua has been termed the "balance 
requirement" by Davis (1963), and represents the excess 
amount to be disposed of from a latitude belt by poleward 
flux-divergence of heat. Table XV(e) shows that the global 
annual average of this parameter (as calculated in this 
study) is zero. 

Warming of the atmosphere by the absorption of solar 
radiation is less than the long-wave cooling rate. The net 
result is that m atmosphere is cooling radiatively at all 
audes ata rate of about 1.13906 55: The compensating 
energy requirement for possible balance has been indicated 
by Equation(76)as arising by sensible and latent heat trans- 
fers across the earth surface. 

With the addition of sensible and latent heat into each 
latitude band, Table XI(e) shows that latitude bands equa- 


torward of 40N have a total energy surplus while latitude 
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bands poleward of 40N exhibit an energy deficit. While the 
mechanisms of horizontal transport are not specified in 
this study, the resultant globally-weighted values of QUA 
of Table XV(e) indicate a mean annual total heat balance 

in the atmosphere. The poleward redistribution required to 
effect this total heat balance is evidently brought about 
by general circulation processes, which were not considered 
in this study. For our purposes it suffices to show that 
the energy-generating aspects of the model proposed here is 


subject to an atmospheric balance. 
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TA, SURFACE HEAT TRANSFERS 


The earth's surface heat balance was computed using 
evaporation and sensible heat transfers over the oceans 
as surface heat sinks. Heat transfer due to evaporation 
and sensible heat are more easily computed over ocean 
surfaces using climatological data available in the Paci- 
fic Ocean representative of the latitude belt as a whole., 
In effect, London's atmospheric soundings were used in 
association with data for ocean surface temperatures and 
surface winds taken from U. S. Weather Bureau's Atlas of 
mn Lic Charts of the Oceans (1938). Continuity of 
sounding data with the oceanic reports was maintained by 
maintaining the same temperature lapse rate in the layer 
zero to 1 km. 

The surface radiative budget is shown schematically 


below: 


Er 7777057 | (77) 


Here Quo is the required oceanic heat transport divergence. 
The surface layer evaporation heat transfer rate was calcu- 
lated after Langlois and Kwok (1969) using: 


2 i = 78 
DT on ln Merle a) (78) 
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277020 the constants: 


2001 5 


Dun = 155 .... gram m 


L 


Also Vo is the climatological ocean surface wind speed in- 
cremented by 2.2 mps in order to account for the effect 

of gustiness upon evaporation [Langlois and Kvok (1969)]. 
A variable drag coefficient Cn proposed by Deacon and Webb 


(1962) for the case of neutral stability was employed: 


1. 4073 
Cn = (1.0 + 0,07 Vg? x 10 


(79) 
At the ocean surface, the atmosphere was assumed to 

be saturated with respect to a plane water surface. The 

surface specific humidity, qun was computed using the 


Clausius-Clapeyron equation: 


622 o AE-BE/TG) 


Ise E me (AE-BE/TG) (80) 


Me specific humidity at the top of the constant flux sur- 
5 — layer, q. ° is given following the parameterization of 
the planetary boundary layer by Langlois and Kwok, by the 


formula: 
HEMOS 0 ТОМ (81) 
0 D 
with: 
A n (82) 
a (95 - вс) 
ии, 
* * 
ИТИ), k and a are turbulent transfer coefficients 


whose values have been documented by Kaitala (1972). 
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li copPree oceans exist equatorvard of 60N for all seasons. 
Therefore, air and sea temperature and wind data were re- 
read for only the first six latitude bands. Since Lon- 
don's atmospheric data was distributed at one-kilometer 


intervals, (Zo ) was taken as 1 km; also T, and qg are 


“10 9 
temperature and specific humidity at the one-km height. The 
first kilometer temperature lapse rate was computed using 
Bondon's soundings. Air temperature at the ocean surface 
was assumed equal to the ocean surface temperature. The 
lapse rate was then applied to yield a nev To for use in 
(82). The new specific humidity qo at 1 km was generated 
using the relative humidity from London's data and the 
Clausius Clapeyron equation. 

Sensible heat transfer rates equatorward of 60N were com- 
puted using Bowen ratio (r) values from Sellers (1965) and 
the formulation: 


Qo = r Qk. (83) 


Poleward of 60N, combined evaporation and sensible heating 
rates were taken from an assessment of meny thousands of 
ship (and also of ice island) reports made by Vowinckel and 
Taylor (1964). 

Tables XVI(a-d) contain the seasonal surface radiative 


balances, which is the oceanic heat transfer require- 


Qvo 3 
ment by poleward flux-divergence. Seasonal variations can 


be compared against the annual mean values listed in Table 


XVI(e). 
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The physical meaning of a negative Qvo value for a 
Etude belt in the tropics (with positive values in 
northerly belts) is that a flux-convergence (divergence) 


of sensible heat due to ocean currents is needed to 


maintain a long-term balance. 
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X. CONCLUSION 


With the use of mean seasonal atmospheric data, a pro- 
posed radiation-heating model for the Naval Postgraduate 
School primitive-equation predictive system was evaluated. 
The model proved readily adaptable to computerization. 

The time and space averaged radiational quantities 
showed rather good agreement as a whole when compared with 
the results of other researchers. Minor difference in the 
comparisons can generally be ascribed to the treatment of 
global albedo, and the assumption of but a single cloud 
layer. While the average global albedo determined in this 
study is typical of usual values quoted, namely 34%, the 
meridional distribution of global albedo was more nearly 
uniform with respect to latitude than would have resulted 
if the mean value had been based upon hourly computations 
for the months involved. This result led to a lower than 
normal value of solar radiation absorbed within the earth's 
surface in low feet dee and a larger than typical value 
absorbed in high latitudes. 

Of the three heat budgets discussed, results compatible 
with Houghton (1954) and Davis (1963) respectively, were 
demonstrated for the earth-atmosphere system and for the 
Meni tude and direction of the required total heat flow. 
Similar tests were also applied to the ocean surface 


where the distribution of surface heat-flux divergence was 
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Enposrte to that expected (e.g. Budyko (1955)), the results 
of the surface heat budget were somewhat biased by the 
averaging involved in the derivation of the global albedo, 
as explained previously. 

Finally, while the radiation package has proven general- 
ly adaptable to a general circulation predictive model, fur- 
ther refinements are anticipated in the treatment of cloud 
cover. A two-layer cloud model for the treatment of middle 
and low-level clouds has been suggested. These two cloud 
types would have different short-wave flux reflection and 
absorption characteristics and would lead to a more accurate 


determination of global albedo. 
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APPENDIX A 


(P Qus TD of the final temperature-dependent flux integral 


€ (u; T,) dB for a variety of final temperature values, (T.) 
and final scaled water-vapor values (ly day”! 


) 


Scaled 
Water 


O 


EF 


195” 


2009 


O 


205 


2109 


215: 


O 


220 


Vapor 


99 


0.0001 26.98 26.53 5051: 51:70 33.45 77. 50077 
0.0005 42.62 45.36 48.19 si, 54.16 31.30 60.55 
0.001 50.84 54.25 57750 61.48 65.30 69.26 7555 
0.002 2.59 64.09 68.46 73.08 27705 82.66 62776 
0.004 59.65 74.74 80.05 85.60 21: 97.41 103.68 
0.006 25.63 81.29 87.20 93239 09.85 106759 110201 
0.008 5:97 86.05 92.41 99.08 106.02 7113.32 12051 
0.01 59737 577 96.52 20527 110.94 118.64 126.69 
0.015 29.63 960.68 7104.09 111.86 120200, E209 55 137.44 
0.02 94.09 101.61 109.52 Jy oS 126.54 125707 145.22 
0.03 100.38 108,57 117,215 9012630 135.83 145.84 156.34 
0.04 104,80 113.49 17705 132.30 142.45 Дези 164.29 
0.06 1210.94 120.33 130.25 140.71 151.73 J 1632 1753351 
0.08 (55.20 125209 155.50 1465000815522 5 17707232 183.42 
2l 118.42 12571 100 15]. 1081608225 176.04 1:5: 
0215 197.08 135.08 146.73 159.07 72217 125,87 200797 
0.2 827.90 139.40 151.59 164.52 178.]9 92:64 20767 
0.3 12.96 145,15 158.10 171.84 186.40 201.81 213; 10 
0.4 136.30 1 162.43 17657 17: :20:206.00007 27.01 
0.5 138.72 151.74 165.61 SO SOS 21» 160982950 F6 
0.6 140.59 105 168.07 183.16 199.19 216.20 234.22 
0.8 143.31 1537.05 m3. 187.34 203.95 221757 240.30 
"о 145.22 159729 174.32 190.33 1 7290 58  225:507, 7244.73 
1.0 148.23 162.86 178,50 195.20 8 1 7. 00752317296 252210 
50 149.99 164.97 11:07 10227 217775 7676 
7: Eu 166.367 182.70 200.18 R212795 T238./6 259507 
3.0 151.89 ASS 183.90 20170 00660 07 07077 70 07 6 SS 
4.0 192.02 168.55 185.44 203355 92222925 9243-61 265.68 
5.0 15.29 100223, 2186.35 204,72 0800 39 245.41 267296 
6.0 153.53 109261 186.90 205.46 225439 246.63 269.40 
8.0 153.99 10088 52157.41 206.24 226.45 218.097 271.23 
EO. O 153.40 169781 1972000 206751 22657: 2021110 0775 
15.0 152.49 LOSOS rT 186,98 206.26 2777007 9. 26 270711 
5.0 150235 167.04 185.09 204 00.1 2/2.26 
30.0 149.32 166.02 184.09 203.61 224.66 247.30 271.63 





APPENDIX B 


Peer vapor equivalent u_ for droplet absorption in relation 
to the pressure-scaled S water-vapor absorber mass u(4,8) 
u le sounding layer (4,8)). Case u (1) corresponds to 
ib 7 titude vater-drop absorptivityş “u (2) to mid-latitude 
absorptivity; u (3) to high latitude absorptivity. 


um. 8) u. (1) u (2): u (3) 

0.01 (1102 102 0 0.05118 0.00740 
0.02 0.18936 0.97292 0.01155 
0.05 0.28886 0.12003 OS 
0.10 0.40852 0.17851 0.03329 
0.15 0.50529 0.22673 0.04365 
0.20 0.56995 0.26939 0.05294 
0.30 0.73785 0.34473 0.96955 
0.50 0.98617 0.47287 0.09826 
0.80 1.29740 0.63548 0.13520 
1.00 1.48119 0.73226 0.15739 
1.20 1.65219 0.8221 0.17822 
1.40 3561292 0.90822 0.19800 
1.60 1.96648 0.98974 0.21691 
1.80 DAD Se 1.06793 0. 23510 
2.00 2.25403 1.14331 DE 
2.20 2.39016 1.213 0.26970 
2.40 252206 1.28699 0.28626 
2.60 2.65022 1.35584 0.30240 
50 2.77499 170207 0.31816 
3.00 2.89675 1.48854 0.33357 
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APPENDIX D 
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The distribution of average temperature and pressure (from 
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